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ABSTRACT 

We present moderate-resolution Keck spectra of nine candidate globular clusters in 
the possible merger-remnant elliptical galaxy NGC 3610. Eight of the objects appear 
to be bona fide globular clusters of NGC 3610. We find that two of the clusters belong 
to an old metal-poor population, five to an old metal-rich population, and only one to 
an intermediate-age metal-rich population. The estimated age of the intermediate-age 
cluster is 1-5 Gyr, which is in agreement with earlier estimates of the merger age, and 
suggests that this cluster was formed during the merger. However, the presence of five 
old metal-rich clusters indicates that a substantial number of the metal-rich clusters in 
NGC 3610 likely came from the progenitor galaxies, although the global ratio of old to 
intermediate-age metal-rich clusters remains very uncertain. 


1. Introduction 

NGC 3610 has attracted attention for some time because of its rich fine structure and promi¬ 
nently warped disk, which imply that a significant dynamical event occurred during the past 5- 
6 Gyr. With a much longer time-scale, relaxation would have led to the disappearance of the fine 
structure in the inner parts of the galaxy. Both the enhanced line strength of H/3 (Schweizer et al. 
1990) and the bluer-than-average UBV colors (Schweizer & Seitzer 1992, hereafter SS92) suggest 
that NGC 3610 harbors an aging starburst population from several Gyr ago, yet the integrated 
colors of the galaxy are not blue enough for it to have undergone a significant burst of star forma¬ 
tion during the past 1 Gyr or so (Scorza &: Bender 1990). NGC 3610 is thus an excellent candidate 
for an intermediate-age merger remnant, with an estimated age of 4 zb 2.5 Gyr (SS92; Whitmore 
et al. 1997, hereafter W97). There is increasing evidence that major bursts of star formation are 
typical among gravitationally interacting galaxies, a good recent example being The Antennae 
(Whitmore et al. 1999). 
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SS92 suggested that elliptical galaxies with intermediate-age stellar populations could be ex¬ 
plained by the occurrence of a major disk-disk merger in their past, which would lend support to 
the theory that some elliptical galaxies formed in such major mergers (Toomre &: Toomre 1972; 
Schweizer 1986; Ashman & Zepf 1992; Barnes 1998). The recent slew of evidence for the formation 
of young massive clusters (YMCs) in many interacting galaxies (see Whitmore 2002 and Schweizer 
2002 for recent reviews) has led naturally to speculations that YMCs are younger versions of the old 
globular clusters seen in the Milky Way and other galaxies. If so, a population of intermediate-age 
clusters could shed light on the suspected merger and might represent the “missing link” between 
YMCs and present-day old globular clusters (GCs). Few systems of intermediate-age GCs are 
presently known—probably the best example being that observed in the merger remnant NGC 
1316 (Goudfrooij et al. 2001a, b). Using a combination of optical/IR photometric data and spectra 
of its globular clusters, these authors were able to estimate, with unprecedented precision, a merger 
age of 3.0 ± 0.5 Gyr. 

W97 reported the discovery in NGG 3610 of a population of GGs that was brighter, redder, 
and more centrally located than the old, metal-poor population. They speculated that this new 
subpopulation might represent a set of clusters formed during a major merger 4 ± 2.5 Gyr ago. In 
order to test this speculation, in the present paper we analyze spectra of nine GGs in NGC 3610 
and determine abundances, metallicities, and ages. 


2. Observations and Data Reduction 

The nine GC candidates were observed with the Keck I and II telescopes and the Low- 
Resolution Imaging Spectrometer (Oke et al. 1995) in multislit mode. Spectra were obtained 
during two observing runs on 1999 April 7-8 (Keck II) and 2000 February 27-28 (Keck I), respec¬ 
tively. Seven target objects for the 1999 run were chosen from the list of GC candidates identified 
by W97 on images taken with the Wide Field/Planetary Camera 2 on the Hubble Space Telescope 
(HST/WFPC2). Two additional objects from that list were added for the 2000 run, based in 
part upon inspection of new, dithered WFPC2 observations described in Whitmore et al. (2002, 
hereafter W02). 

During the 1999 run, ten exposures totaling 300 minutes were taken with a single mask over 
two nights. During the 2000 run, nine exposures totaling 260 minutes were taken with a new mask, 
also over two nights. All observations were made with a 600 line mm“^ grating blazed at 5000 
A. This grating provided a reciprocal dispersion of 1.28 A pixeU^ and a spectral resolution of ~6 
A. The total wavelength range covered by the CCD detector was generally 3600 — 6000 A, but 
the usable range, especially at the blue end of the spectrum, was often shorter by several hundred 
Angstroms due to the differences in the location of individual slitlets within the mask. 
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The standard data reduction was performed using IRAF^ and the data analysis package BAO- 
LAB, written by one of us (SL). Raw images were debiased and then flatfielded, using a normalized 
composite sky flat. After cosmic-ray removal using the IRAF task cosmicrays, the BAOBAB tasks 
xdistmap and xdistapp were used to rectify the science spectra, since several of the slitlets were 
slightly tilted with respect to the spatial axis. This tilting was done to allow more cluster candi¬ 
dates to fit on the slitmask. The extracted spectra were then wavelength calibrated using spectra 
of comparison arc lamps. 

After this wavelength calibration, the science spectra were shifted by small amounts in wave¬ 
length in order to bring several bright sky lines, principally those of oxygen, into agreement with 
their laboratory wavelengths. The spectra were then coadded with sigma clipping. The coadded 
spectra were flux-calibrated using observations of the flux standard PG 0823-1-546, chosen from 
Massey et al. (1988). To assess the quality of the observations, the signal-to-noise ratio (S/N) of 
the spectra was measured in the range 4700 — 5300 A, where most features of interest lie. 

Table 1 lists basic data for the nine clusters in our sample. The positional and photometric data 
are taken from W02, while the heliocentric radial velocities are based on our own measurements 
described in §3. No extinction corrections were applied since Burstein & Heiles (1984) give Ay = 
0.00 mag and Schlegel et al. (1998) give Ay = 0.03 mag. 

One item of note is that during the 2000 run, the important H/3 feature in the spectrum of 
object W3 was marred by a block of bad pixels. Since no effective correction could be applied, 
our subsequent analysis of W3 includes data solely from the 1999 run. This problem is not overly 
disconcerting, however, since the quality of the 1999 spectrum alone is relatively high (S/N ~ 25). 

Spectra of a representative sample of clusters (W6, W9, and WIO) are presented in Figure 1. 


3. Radial Velocities 
3.1. Velocity Measurements 

The radial velocities of individual clusters were measured by cross-correlating the cluster spec¬ 
tra with spectra of the radial velocity standards HR 1805 and HR 3905 using the IRAF task fxcor. 

For cluster W26* no clear cross-correlation peak with the standard stars could be found. While 
a reasonable value was obtained by cross-correlating the cluster spectrum with spectra of several 
other clusters, the correlation coefficient of the match was low, indicating the measured radial 
velocity may have been spurious. Though these difficulties could be due to the relatively low S/N 
of the cluster, comparison of the smoothed spectrum of W26* with those of W6 and WIO revealed 
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little resemblance. The spectrum of W26* shows none of the strong lines (e.g. H/3, Mg) present in 
the other clusters, nor does it show an H+K line break. Furthermore, the strength of its spectrum 
increases towards the UV, unlike those of the other clusters. Visual examination of the WFPC2 
images show that the cluster appears superimposed on an extended red object, which may be a 
background galaxy. If indeed the cluster is real, the observed spectrum may be dominated by that 
of the galaxy. For these reasons, we have excluded W26* from further analysis, and it does not 
appear in any of the subsequent tables or figures. 

The mean heliocentric radial velocity of the eight clusters is 1769 ± 15 km s“^, with a velocity 
dispersion of 45 ± 25 km s“^. This mean velocity is in good agreement with the value 1787 ± 48 km 
s“^ for the systemic velocity of NGC 3610 from the Third Reference Catalogue of Bright Galaxies 
(de Vaucouleurs et al. 1993), but does differ significantly from the more recent value of 1696 ± 17 
km s“^ listed in the Updated Zwicky Catalogue (Falco et al. 1999). 

There is no clear radial trend in the eight cluster velocities, nor do there appear to be trends 
with V magnitude, color, or (based upon data in §4) [Fe/H]. There are no distinct kinematic 
subgroups visible either. With such a small sample size, these null results are not unexpected. 


3.2. Galactic Mass 

The estimated errors of the radial velocities are small enough that we can make a rough estimate 
of the mass of NGC 3610 lying within the orbits of our cluster sample. We use the new tracer mass 
estimator (Evans et al. 2002), which has numerous advantages over the commonly-used projected 
mass estimator (Bahcall &: Tremaine 1981; Heisler et al. 1985). These include; no assumption that 
the tracers follow the dark matter distribution, the use of an inner and outer radius for the tracer 
population, and a smaller dependence on the exact orbital shapes. (Evans et al. find that assuming 
isotropic orbits only introduces a 30% uncertainty in the final mass estimate). If we assume that 
the potential is isothermal and the GC population is isotropic and falls off as ~ r~^, then the 
equation for the tracer mass estimator is; 


M = 


16 (rin/rput) ^ - 1 

ttGN log(rout/nn) 


N 




( 1 ) 


where and Tout are the inner and outer radii, respectively, of our tracer population. We find that 
M = 8.6 X 10^^ Mq within a radius of 13 kpc. The sources of error in this estimate are numerous, 
including measuring errors in the radial velocities, the limited sample size, our lack of knowledge of 
the true GC orbits and radial distribution, and the assumption of an isothermal potential (especially 
in the context of the inner regions of a merger remnant). The uncertainties are great enough that 
it makes little sense to give a formal error estimate, whence we will express our mass estimate as 
8 X 10^*^ Mq, with the caveat that this value is uncertain by at least a factor of two. 


r\j 
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4. Metallicities, Abundances, and Ages 
4.1. Metallicities 

Brodie & Huchra (1990, hereafter BH90) defined a procedure for measuring the metallicity of 
globular clusters by taking a weighted mean of various elemental absorption-line indices sensitive 
to [Fe/H], Indices defined over a wavelength range of interest are calculated with respect to a 
pseudocontinuum, defined using regions to either side of the feature passband. 

Before measurement of any features, the spectra were flux calibrated as described in §2. The 
minor role of flux calibration in measuring indices has been documented extensively (Faber et ah 
1985; Kissler-Patig et ah 1998; Larsen & Brodie 2002). Errors due to flux calibration are extremely 
small, amounting to ~5% of the smallest error bars in our study. No resolution corrections were 
made for measuring BH90 metallicity indices, since the BH90 calibration primarily used data near 
the resolution of our data or better. The spectra were then corrected for radial-velocity shifts using 
the values given in Table 1. 

We measured [Fe/H] for our eight clusters according to the formulae in BH90, with the excep¬ 
tion of a slight modification in calculating the variance, as described in Larsen & Brodie (2002). 
The resulting individual index [Fe/H] estimates are listed in Table 2. 

To provide a check on these spectroscopic metallicity estimates, we also calculated photometric 
[Fe/H] values according to the Kissler-Patig et al. (1998) equation: 

[Fe/H] = (-4.50 ± 0.30) + (3.27 ± 0.32)(K - I) (2) 

The photometric and spectroscopic [Fe/H] values are given in Table 3, and the estimates 
compared in Figure 2. Seven of the eight clusters fall near a 1:1 relation, indicating good agreement 
between the two value. The one outlier is W30, which has a low-S/N spectrum. The supersolar 
spectroscopic [Fe/H] value for this cluster is particularly suspect because of the wide spread in 
the individual metallicity indices, which range from [Fe/H]^ = —2.00 to [Fe/H]^^^^ = 3.66. This 
suggests that while it may be possible that the cluster has a metallicity near solar (such clusters 
have recently been discovered in NGC 1399 by Forbes et al. 2001), the large error bars make it 
impossible to put any strong constraints on its metallicity at this time. 

Setting aside the cluster W30 due to its large error bars, and using a dividing line of [Fe/H] = 
— 1.00 between metal-rich and metal-poor clusters, the mean metallicities of these two subpopula¬ 
tions are [Fe/H] = —0.69T0.11 and —1.28T0.15, respectively. These values are certainly consistent 
with the “standard” values of —0.5 and —1.5 (Larsen et al. 2001; Kundu & Whitmore 2001), espe¬ 
cially considering the small-number statistics operating. There is generally a great deal of variation 
in the location of the metal-rich peak among galaxies, and Larsen et al. (2001) find that the scatter 
in the color of the metal-rich peak is 0.04 mag, equivalent to ~0.13 dex in metallicity, using the 
Kissler-Patig et al. (1998) relation. 
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4.2. Lick/IDS Indices 

Lick/IDS indices have evolved from the original work of Burstein et al. (1984, hereafter B84), 
with several authors providing refinements to the system, usually in the form of additional indices 
that cover new features or define the pseudocontinuum regions of an existing feature passband dif¬ 
ferently (Worthey et al. 1994; Trager et al. 1998). These refinements have been accompanied by new 
sets of cluster-evolution models that predict how various Lick/IDS indices (e.g., H/3, Mg 2 , Fe5335) 
evolve in age-metallicity space (Bruzual & Chariot 2002; Maraston & Thomas 2000; Maraston et 
al. 2001; Worthey et al. 1994; Alvensleben & Burkert 1995). The refined spectroscopic indices 
and new models are especially important in the study of intermediate-age clusters, since the colors 
of a 3-4 Gyr old population can easily blend with those of a much older metal-rich component, 
rendering photometric data insufficient to determine ages (Worthey et al. 1994; W97). In passing 
we note that this degeneracy can be partially lifted through the use of infrared photometry. For 
example, Puzia et al. (2002) demonstrate that subpopulations which separate cleanly in a VIK 
color-color diagram would blend together in the common VI color-magnitude diagram. 

For the measurement of Lick/IDS indices we had to correct for spectral resolution. Since the 
calibration data came from sources with a spectral resolution range of 5-12 A, we convolved our 
spectra with a 5-pixel (equivalent to 6.4 A) Gaussian kernel as a “middle-of-the-road” approach to 
more closely match the resolution of the original Lick data. Although some previous studies (e.g., 
Kissler-Patig et al. 1998) have found little or no resolution degradation necessary when compared 
to the standard values in Worthey et al. (1994), other studies (e.g., Larsen & Brodie 2002) have 
found offsets of 0.2 — 0.3 A with respect to those standard values when measuring indices in 
unsmoothed LRIS spectra. These differences are comparable to the error bars on many of our 
index measurements. Thus, we decided that smoothing our spectra was the best course of action. 
It is, however, worthwhile to recall the Kissler-Patig et al. (1998) warning that “care should be 
taken in the intercomparison of measurements of various groups.” 

For this study, we measured a set of Lick/IDS indices for comparison to theoretical models. 
These included the higher-order Balmer indices and H5a (Worthey & Ottaviani 1997). Table 
4 gives these indices for all eight clusters. 


4.3. Ages 

To estimate the ages of the GGs in our sample, we plotted various combinations of the age- 
sensitive indices (e.g., H/3 and Hy) against the metallicity features (e.g., Mg 2 , Mgb, Fe5270, and 
Fe5335). As an example. Figure 3 shows H/3 plotted against the averaged Fe5270-|-Fe5335 index. 
Superimposed on this plot are Maraston theoretical isochrones and isometallicity lines. 

Figure 3 clearly shows three distinct subgroups of clusters; an old metal-poor population, an 
old metal-rich population, and a single metal-rich cluster (W6) of intermediate age. All of the old 
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clusters fall below the 14 Gyr isochrone, but the model age uncertainties are great enough that this 
is not particularly worrisome. Indeed, a general feature of simple stellar population (SSP) models 
is that while the absolute ages are uncertain, relative ages are usually consistent. Within the error 
bars, all seven old clusters seem to be coeval. 

While the Maraston (and other SSP) models are intended mainly to represent old stellar 
populations, and thus may be less accurate for intermediate-age clusters, it is clear that the cluster 
W6 cannot be older than ~5 Gyr. The absence of strong Balmer lines implies that the cluster also 
cannot be < 1 Gyr old. Therefore, the age of W6 seems reasonably constrained to lie in the range 
1-5 Gyr. A literal interpretation of Figure 3 suggests an age of 2.5lQg Gyr. 

These arguments are confirmed by an H/3 vs. Mg 2 diagram (Fig. 4) and an H/3 vs. [MgFe] 
diagram (Fig. 5). Though the old metal-poor and metal-rich subpopulations appear blurred, the 
general result is the same: one of the clusters is relatively young, while the rest are old and—within 
the error bars—coeval. (No confirmation could be made using the age-sensitive indices and 
H5y4, as their error bars are too large.) One advantage of using the combined [MgFe] index is that 
it is insensitive to variations in [a/Fej. 

The absolute magnitude of cluster W6 is consistent with it being relatively young; at My = 
— 10.60 (W02), it is the sixth-brightest cluster in NGC 3610. This is more luminous than any Milky 
Way GG, and nearly as luminous as the brightest GC in M31, Gl, which has My = —10.85 (Reed et 
al. 1994) (However, Barmby et al. (2002) have recently presented evidence that the heavily reddened 
cluster 037-B327 may have My = —11.74). Excepting the outlying [Fe/H] value of cluster W30, 
W6 is also the most metal-rich GG in our sample with [Fe/H] = —0.29 ±0.19. Note, however, that 
the BH90 calibration is based on a sample of old galactic and M31 GCs, whence this metallicity 
estimate could be significantly less accurate than the formal error bars indicate. Lick/IDS Fe and 
Mg indices for W6, interpreted with Maraston models (Figs. 3 &: 4), indicate the metallicity of the 
cluster may be solar or supersolar, much higher than the previous estimate. 

Our age estimate for cluster W6 is in good agreement with the merger age of 4±2.5 Gyr deduced 
by SS92 and W97, and provides evidence in favor of the idea that W6 is an intermediate-age object 
formed in a merger several Gyr ago. 

Figure 6 shows a color-magnitude diagram of all GCs detected in the W02 survey, with our 
eight clusters (and the intermediate-age cluster W6) specially marked. 


4.4. [a/Fe] Abundances 

To estimate [a/Fe] for the clusters in our sample, we used the new a-enhanced models of 
Thomas et al. (2002). Figure 7 is an Mgb vs. <Fe> diagram for the seven old GCs, overplotted 
with 14 Gyr model isochrones for [a/Fe] = 0.0 to ±0.5. Though there is substantial scatter among 
the points, there seems to be a hint of a spread in [a/Fe], ranging from ~ 0.0 to ±0.3. The solar a- 



element abundance ratios for several old clusters suggest that either the GCs formed over a period of 
several Gyr, during which time there was substantial enrichment of the ISM by Type la supernovae 
(SNe), or the GCs with differing a-enhancement came from different progenitor galaxies; within 
each individual galaxy, GC formation took place on a short timescale. The latter scenario offers 
the intriguing possibility of associating individual GCs with their parent galaxies. 

Figure 8 shows a similar diagram for the intermediate-age GC W6, with 3 Gyr isochrones 
overplotted. Interestingly, W6 seems to have [a/Fe] ~ -1-0.3. A merger-induced starburst, created 
as gas from each of the progenitor disks is funneled to the center of the potential well, would 
be expected to enrich the ISM with a-elements through Type II SNe within a few hundred Myr. 
Assuming that the gas started with [a/Fe] ~ -1-0.0, this enrichment would have been substantial, 
raising the absolute a-element abundance of the ISM from half-solar to solar. If so, W6 must then 
have formed before the Type la SNe had time to bring [a/Fe] back up to solar, suggesting that the 
cluster was formed in the early days of the merger. 

Alternatively, if W6 formed with a very top-heavy IMF, internal enrichment could be respon¬ 
sible for the supersolar [a/Fe] ratio. Whether such GCs exist in other galaxies is unclear; Brodie 
et al. (1998) and Larsen & Brodie (2002) suggest that the anomalously strong Balmer lines in a 
few young GCs in NGC 1275 and NGC 1023, respectively, could be caused by a top-heavy IMF. 
However, the significance of these results is colored by the uncertainties in SSP models at such 
young ages. 


5. The Nature of the Merger 

Even before GCs were studied in NGC 3610, several pieces of evidence pointed toward this 
E5 galaxy being the remnant of a major merger. Both its unusually rich fine structure (Seitzer 
& Schweizer 1990) and warped central disk (W97) strongly suggest that it is dynamically young. 
Given the estimated merger age of ~4 Gyr, the presence of sharp-edged ripples (“shells”), luminous 
plumes, and an exceptionally boxy halo (Schweizer 1998, esp. Eig. 41) can hardly be attributed 
to the infall of a minor companion many orbital periods ago. Rather, such a halo shape and 
dynamically cold features tend to be signatures of merged disk galaxies and of the delayed return 
of their tidally ejected material (see Barnes 1998 for a review). 

Supporting the morphological evidence in this otherwise normal small-group elliptical is strong 
evidence for an aging starburst; spectral indices indicate enhanced H/3 absorption and weakened 
Mg and CN features, as expected from an aging starburst superposed on old stellar populations 
(Schweizer et al. 1990). From the bluer-than-average UBV colors and a heuristic merger model, 
SS92 estimate a likely merger age of around 4 Gyr. Near-infrared JHK photometry and HST/STIS 
spectra in the space UV indicate the presence of a significant intermediate-age population as well 
(Silva & Bothun 1998; Spinrad 1997). Gontrary to the claims by Silva & Bothun, the presence of 
a central stellar disk is perfectly compatible with a major disk-disk merger origin of this galaxy; 
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gas-rich mergers tend to form first gaseous, and then stellar, central disks from tidally ejected 
and returning gas (Hernquist & Barnes 1991; Barnes 2002), as directly observed in some nearby 
remnants of recent mergers. 

The study of GCs in NGC 3610 adds interesting new information plus some question marks 
to this proposed disk-merger formation scenario. 

Like the majority of elliptical galaxies, NGG 3610 harbors a GG system with a bimodal color 
distribution indicative of at least two cluster subpopulations (W97; W02). The spatial distribution 
of the red subpopulation appears more centrally concentrated than that of the blue subpopulation, 
suggesting that additional gaseous dissipation occurred between the formation of the metal-poor 
and metal-rich clusters. What is unusual about the GC system of NGC 3610 is that the brightest red 
GCs are about 0.7 mag more luminous in V than the brightest blue clusters (though spectroscopy 
is necessary for certain discrimination of metal-rich and metal-poor clusters). This magnitude 
difference, together with the measured color difference A{V — I) between the blue and red peaks 
of the color distribution, suggests that at least some of the red GCs may be only about 4± 2.5 Gyr 
old (see Figs. 13 &: 15 in W97). 

The luminosity functions of the two GC subpopulations appear to support a significant mean- 
age difference: while the luminosity function of the blue GCs is approximately lognormal, that of 
the red GCs is well represented by a power law with index a ~ —1.8 (W02). This power-law shape 
suggests that the red subpopulation is dynamically much less evolved and, hence, significantly 
younger than the blue subpopulation. (However, see below for a caveat on the red LF). 

Our spectroscopic observations now add more detailed age and metallicity information, though 
only for eight GCs, all lying beyond a projected radius of 3 kpc from the center. 

First, one of the two sampled GCs nearest the center appears indeed to be of intermediate age, 
with an estimated age of 1-5 Gyr. This strengthens the cases previously made for the presence of 
stellar and GC populations of intermediate age in NGC 3610. According to the BH90 calibration 
(Table 3) and the Lick/IDS indices interpreted with Maraston models (Figs. 3 & 4) the metallicity 
of this one cluster is [Fe/H] > —0.3 (see §4.3), while the Mg abundance may be solar or slightly 
higher (Fig. 4). The metallicity and Mg abundance argue strongly against this cluster having 
formed from the metal-poor gas of an infalling dwarf galaxy. For example, the mean metallicity of 
14 LMC clusters in the age range 1.5-4 Gyr is only ([Fe/H]) = —0.74, with a scatter of o" = 0.22 
(Leonardi 2000). Hence, any LMC-like dwarf falling into NGC 3610 1.5-4 Gyr ago would seem 
unlikely to contribute or produce clusters of the metallicity observed in W6. Rather, a metallicity 
of [Fe/H]> -0.3 is in accord with the gaseous metallicities that one might expect in some average 
Sb or Sc spirals about 2-4 Gyr ago (e.g., Pagel & Tautvaisiene 1995). 

Second, both Figs. 3 and 4 suggest that there is a significant age difference between the one 
intermediate-age cluster W6 and the other five “metal-rich” clusters, the latter apparently being 
older than ~10 Gyr. In accord with this age difference, the metallicities of these five clusters are 
significantly lower as well, averaging perhaps around [Fe/H] ss —0.7 to —1.0 (see also Table 3). 
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Thus, the red subpopulation of GCs in NGC 3610 seems to contain a broad range of cluster ages 
and metallicities. Specifically, the ages of the six red clusters suggest that there were at least two 
major formation episodes of metal-enriched GCs. This raises a currently unanswerable question: 
What fraction of the red GCs formed during the suspected major-merger event ~4 Gyr ago? This 
is an important factor in determining how the specific GC frequency in this galaxy will evolve over 
a Hubble time. 

With our present sample of only six red clusters, of which one appears to be of intermediate 
age, the relative numbers of old and intermediate-age objects are poorly constrained. Furthermore, 
the spectroscopic sample is clearly biased toward the outer halo, while the intermediate-age GCs 
may be preferentially located in the inner parts of the galaxy. 

While the power-law luminosity function of the red GCs would seem to indicate that many 
of these clusters may be of intermediate age rather than old, one would clearly expect a number 
of metal-rich, old clusters from the progenitor galaxies. In the Local Group (dominated by the 
GC systems of the Milky Way and M31), the overall ratio of metal-poor to metal-rich GCs is 
about 2.5:1 (Forbes et al. 2000). Thus, we might expect at least 30% of the old GCs in NGC 
3610 to be metal-rich, although the number of metal-rich clusters probably depends on the Hubble 
type of the progenitor spirals. Comparison of the relative numbers of old and intermediate-age 
objects is further complicated by their different luminosity functions, age fading, and the fact that 
a portion of the intermediate-age objects are going to disintegrate as their LF evolves towards a 
more nearly log-normal shape (Fall & Zhang 2001; Vesperini 2001). Spectroscopy of a larger sample 
of GCs extending to smaller radii, where—unfortunately—the galaxy background becomes a major 
obstacle, will be necessary to further constrain the relative numbers of old and intermediate-age 
objects. 


6. Conclusions 

We have observed nine globular cluster candidates in the disturbed elliptical galaxy NGC 3610, 
which is thought to be a merger remnant. Eight of the nine objects appear to be bona fide globular 
clusters belonging to NGC 3610. Seven of these clusters fall into the common old metal-poor and 
old metal-rich subpopulations. However, there is considerable evidence for an intermediate-age 
metal-rich cluster that likely formed in the merger event 1-5 Gyr ago. From the cluster kinematics 
and using the tracer mass estimator, we estimate the mass of the galaxy within a radius of 13 kpc 
to be ~ 8 X IO^'^Mq. Finally, we point out that the ages and metallicities of clusters in the red 
subpopulation suggest at least two cluster formation episodes, although the relative importance of 
each of these two formation epochs in creating the current metal-rich cluster population is unknown. 

We thank Linda Schroder for help with the observations. We also gratefully acknowledge 
support by the National Science Foundation through Grants AST-9900732, AST-0206139 to JPB 
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Fig. 1.— Spectra of clusters W6 (intermediate-age metal-rich), W9 (old metal-rich), and WIO (old 
metal-poor). Notice the strong H/3 line in W6 and the decreasing strength of metal lines from 
W6-W10. 
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. 3.— Plot of H/3 vs. a composite Fe index, with a grid of model isoch: 
IS by Maraston superposed. From left to right, the isometallicity lines re] 
35, -0.33, 0.00, 0.35. 
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Fig. 5.— Plot of H/3 vs. [MgFe], with a grid of model isochrones and isometallicity lines by 
Maraston superposed. [MgFe] is a composite index insensitive to [a/Fe] variations. From left to 
right, the isometallicity lines represent [Fe/H] = —2.25, —1.35, —0.33, 0.00, 0.35. 
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Fig. 6.— A color-magnitude diagram of all GCs detected in the HST survey by W02. The open 
circles are objects detected on the PC chip, and are too close to the galaxy center for spectroscopy. 
The filled circles are objects detected on the three WF chips. Stars mark the positions of the seven 
old GCs in our sample; a triangle marks the intermediate-age GC W6. The dashed line denotes 
the photometric color cut iV — I = 1.025) adopted by W02 to separate blue and red clusters. A 
typical error bar for the clusters in our sample is shown. 
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Fig. 7.— Comparison of Mgb and <Fe> for the eight old GCs. 14 Gyr isochrones from Thomas et 
ah, varying in [a/Fe], are superimposed. 






Table 1. Basic Data for Globular-Cluster Candidates in NGC 3610 


ID ^ 

R.A. 

(hr;min:sec) 

Dec. 

(°:':") 

V 

(mag) 

V-I ^ 
(mag) 

RV 

(km s”-^) 

Proj. Radius 

CO 

'^’‘^Exp. Time 
(min) 

S/N 

W3 

11:18:27.96 

58:47:16.05 

21.53 ±0.05 

1.05 ± 0.05 

1831 ± 20 

22.0 

300 

24.7 

W6 

11:18:28.03 

58:47:05.11 

21.79 ±0.05 

1.24 ±0.05 

1799 ± 17 

22.7 

560 

30.4 

W9 

11:18:33.84 

58:46:53.75 

21.96 ±0.05 

1.15 ±0.05 

1825 ± 22 

69.3 

560 

36.7 

WIO 

11:18:32.95 

58:46:15.41 

22.00 ± 0.05 

0.94 ±0.05 

1774 ± 32 

82.0 

560 

36.5 

W26* 

11:18:33.52 

58:47:42.22 

23.58 ±0.03 

1.22 ±0.04 


72.4 

260 

8.5 

W28 

11:18:24.71 

58:45:53.33 

22.85 ±0.05 

1.07 ±0.05 

1712 ± 32 

77.9 

260 

10.4 

W30 

11:18:30.02 

58:47:28.42 

23.01 ±0.05 

1.09 ±0.05 

1724 ± 31 

41.2 

560 

13.4 

W31 

11:18:26.45 

58:46:29.29 

23.08 ± 0.05 

1.18 ±0.05 

1755 ± 21 

42.9 

560 

12.3 

W32 

11:18:32.75 

58:46:03.85 

23.09 ± 0.05 

1.16 ±0.05 

1732 ± 24 

89.2 

560 

12.9 


“^Taken from Whitmore et al. (2002), except for W26*, which is from Whitmore et al. (1997). 
’^Heliocentric radial velocities determined in this work. 

'^At the distance of NGC 3610, a projected radius of 10" corresponds to 1.46 kpc. 

^Total exposure time over both observing runs. 

^Average signal-to-noise ratio per resolution element over the wavelength range 4700-5300 A. 



Table 2. [Fe/H] Estimates from Individual BH90 Indices 


ID 

A 

Mg2 

MgH 

G Band 

CNB 

Fe5270 

CNR 

H ± K 

W3 

-1.05 ±0.37 

-0.71 ±0.35 

-0.18 ±0.50 

-1.01 ±0.40 

-0.39 ±0.37 

-0.35 ±0.64 

-0.96 ± 0.47 

-1.12 ±0.51 

W6 

-1.16 ±0.37 

0.00 ±0.35 

0.60 ± 0.49 

-0.41 ±0.35 

-0.32 ±0.35 

-0.32 ±0.63 

-0.42 ±0.46 

-.64 ±0.45 

W9 

-1.55 ±0.37 

-0.66 ±0.34 

-0.65 ±0.49 

-1.33 ±0.35 

-1.71 ±0.36 

-0.57 ±0.62 

-1.12 ±0.46 

-1.81 ±0.45 

WIO 

-1.39 ±0.37 

-1.65 ±0.34 

-1.39 ±0.49 

-1.28 ±0.34 

-1.99 ±0.34 

-1.52 ±0.63 

-1.26 ±0.46 

-2.51 ±0.44 

W28 


-1.75 ±0.40 

-0.43 ±0.62 

-0.22 ±0.61 


-0.23 ±0.78 

-2.57 ±0.50 

0.23 ± 0.76 

W30 

-2.00 ±0.37 

-0.53 ±0.39 

0.53 ± 0.58 

0.01 ±0.46 

3.66 ± 0.47 

-0.09 ±0.73 

-0.52 ±0.50 

1.29 ±0.47 

W31 

-2.11 ±0.37 

-1.24 ±0.38 

-1.38 ±0.56 

0.67 ±0.48 

0.98 ±0.52 

-0.49 ±0.77 

-1.02 ±0.49 

-1.85 ±0.51 

W32 

-1.82 ±0.37 

-1.28 ±0.39 

-1.11 ±0.57 

-1.22 ±0.56 

-0.41 ±0.44 

0.12 ±0.72 

-1.30 ±0.49 

-0.12 ±0.50 
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Table 3. Spectroscopic & Photometric [Fe/H] Values 


ID 

[Fe/H] (spec.) 

[Fe/H] (phot.) 

W3 

-0.71 ±0.15 

-1.07 ±0.48 

W6 

-0.29 ±0.19 

-0.45 ±0.52 

W9 

-1.18 ±0.20 

-0.74 ±0.50 

WIO 

-1.65 ±0.18 

-1.43 ±0.46 

W28 

-1.02 ±0.57 

-1.00 ±0.49 

W30 

0.40 ± 0.69 

-0.94 ±0.51 

W31 

-0.77 ±0.48 

-0.64 ±0.50 

W32 

-0.89 ±0.27 

-0.71 ±0.50 




Table 4. Lick/IDS Indices®" 


ID 

H/3 

(A) 

H7a 

(A) 

H5a 

(A) 

Ca4227 

(A) 

G4300 

(A) 

Fe5270 

(A) 

Fe5335 

(A) 

Mg2 

(mag) 

Mgb 

(A) 

CN2 

(mag) 

W3 

1.17 ±0.32 

-3.16 ±0.55 

-0.25 ± 0.73 

-0.23 ±0.38 

3.85 ± 0.60 

2.53 ±0.37 

2.87 ±0.47 

0.16 ±0.01 

2.75 ±0.37 

0.06 ± 0.02 

W6 

2.59 ±0.27 

-2.90 ±0.33 

-1.24 ±0.51 

1.14 ±0.27 

5.16 ±0.36 

2.66 ±0.28 

2.58 ±0.33 

0.23 ±0.01 

3.99 ±0.29 

0.13 ±0.01 

W9 

1.42 ±0.21 

-2.22 ±0.39 

1.33 ±0.38 

0.82 ± 0.22 

2.89 ±0.39 

2.46 ±0.25 

2.22 ±0.30 

0.16 ±0.01 

3.72 ± 0.25 

0.00 ±0.01 

WIO 

2.34 ±0.22 

-0.90 ±0.34 

1.09 ±0.41 

0.24 ±0.19 

3.05 ± 0.34 

0.82 ±0.27 

1.15 ±0.30 

0.06 ±0.01 

0.92 ±0.24 

0.01 ±0.01 

W28 

2.28 ±0.81 

-0.89 ± 1.25 

6.16 ±0.97 

2.61 ±0.70 

5.53 ± 1.21 

3.07 ±0.79 

-2.80 ± 1.12 

0.06 ±0.02 

2.13 ±0.84 

-0.05 ± 0.04 

W30 

1.50 ±0.62 

-3.99 ± 1.06 

0.56 ±0.94 

0.47 ±0.27 

5.43 ± 0.75 

2.57 ±0.67 

2.03 ± 0.84 

0.18 ±0.02 

3.71 ±0.56 

0.11 ±0.03 

W31 

0.92 ±0.63 

-3.82 ±0.90 

5.82 ± 1.09 

-0.18 ±0.73 

5.95 ±0.85 

2.52 ±0.79 

1.62 ±0.97 

0.11 ±0.02 

1.98 ±0.80 

0.01 ±0.02 

W32 

0.93 ±0.62 

-3.14 ± 1.04 

-2.14 ±0.61 

0.55 ± 0.55 

2.29 ± 1.13 

3.22 ±0.64 

1.87 ± 1.03 

0.10 ±0.02 

0.55 ±0.76 

-0.01 ±0.03 


These indices were measured according to the definitions in Trager et al. (1998) and Worthey & Ottaviani (1997). 



